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FAST SYNTHESIS SUB-BAND FILTERING FOR DIGITAL SIGNAL DECODING 

This invention relates to digital signal decoding for the purposes primarily of audio 
reproduction. In particular, the invention relates to enhanced synthesis sub-band filtering 
5 during decoding of digital audio signals. 

In order to store or transmit data representing audio signals it is often desirable to first encode 
or compress the data so as to enable it to be stored or transmitted more efficiently. Decoding 
the data requires that the stored or transmitted data be reconstructed into audio signals by 

1 0 application of a decoding or decompression technique. The reconstruction process is typically 
quite computationally intensive, yet the process should be fast and reliable enough to enable 
the audio signals to be reconstructed in real time, on the fly, for example. In order for the 
decoding process to be carried out in relatively low-cost consumer products, the hardware 
utilised by the decoder should also preferably be relatively simple and mexpensive, or at least 

15 to the greatest extent reasonably possible. 

Efficient stereo and multichannel digital audio signal codmg methods have been developed for 
storage or transmission applications such as Digital Audio Broadcasting (DAB), Integrated 
Service Digital Network (JSDK), High Definition Television (HDTV) and Set Top Box (STB) 
20 for video-on-demand. The formats used to encode and reciprocally decode digital audio and 
video mformation for storage and retrieval is subject to various standards, one of which has 
been established by the Movmg Picttires Experts Group and is known as the MPEG standard. 
A standard on low bit rate codmg for mono or stereo audio signals was established by MPEG- 
1 Audio, published under ISO-IEC/JTCl SC29 11172-3, entitled "Coding of Moving Pictores 
25 and Associated Audio for Digital Storage Media at up to About 1.5 Mbit/s", and the 
disclosure of that document is incorporated herein by reference. MPEG-2 Audio (ISO/IEC 
13818-3) provides the extension to 3/2 multichannel audio and an optional low frequency 
enhancement channel (LFE). The audio part of the standard, ISO/IEC 11172-3, defines three 
algorithms, Layer 1, 2 and 3 for coding PCM audio signals. MPEG-2 (Multichannel) also 
30 defines Layer 1, 2 and 3 algorithms. 
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The MPEG audio encoder processes a digital audio signal and produces a compressed 
bitstream for transmission or storage. The encoder algorithm is not standardised, and may 
use various means for encoding such as estimation of the auditory masking threshold, 
quantisation, and scaling. However, the encoder output must be such that a decoder 
5 conforming to the above-mentioned standards specification will produce audio suitable for the 
intended application. 

The decoder, subject to the application-dependent parameters, accepts the compressed audio 
bitstream in the defined syntax, decodes the data elements and uses the information to produce 

10 digital audio output, also according to the defined standard. The decoder first unpacks the 
received bitstream to recover the encoded audio information frame by firame. After the 
process of frame unpacking, the decoder performs an inverse quantisation (expansion process) 
and feeds a sub-band synthesis filter bank with a set of 32 scaled-up sub-band samples in 
order to reconstruct the output PCM audio signals. The sub-band filter banks used for Layer 

15 1 and Layer 2 of MPEG 1 audio decoder and Layer 1 and Layer 2 of MPEG2 (Multichannel 
extension) audio decoder, are the same. 

The sub-band synthesis filter is one of the most computationally intensive blocks of the MPEG 
audio decoder. Sub-band filtering is performed for each sub-band in a frame and for every 
20 channel. Any reduction in its computational requirements thus enables less complexity and 
reduced cost of decoding. 

In accordance with the present invention there is provided a method of decoding digital audio 
data, comprising the steps of obtaining an input sequence of data elements representing 
25 encoded audio samples, calculating an array of sum data and an array of difference data using 
selected data elements from the input sequence, calculating a first sequence of output values 
using the array of sum data, calculating a second sequence of output values using the array 
of difference data, and forming decoded audio signals from the first and second sequences of 
output data. 

30 

Preferably, the array of sum data is obtained by adding together respective first and second 
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data elements from the input sequence, the first and second data elements being selected from 
mumally exclusive sub-sequences of the input sequence. Furthermore, the array of difference 
data is preferably obtained by subtracting respective first data elements from corresponding 
second data elements of the input sequence, the first and second data elements being selected 
5 from mumally exclusive sub-sequences of the input sequence. 

In one form of the invention the step of calculating an array of sum data and an array of 
difference data comprises dividing the input data sequence into first and second equal sized 
sub-sequences, the &st sub-sequence comprising the high order data elements of die input 

1 0 sequence and the second sub-sequence comprising the low order data elements of the input 
sequence, calculating the array of sum data by adding together each respective data element 
of the first sub-sequence with a respective corresponding data element of the second sub- 
sequence, and calculating the array of difference data by subtracting each respective data 
element of the first sub-sequence from a respective corresponding data element of the second 

15 sub-sequence. 

The invention also provides method of decoding a sequence of m, m an even positive integer, 
input digital audio data samples S[k], where k = 0, 1, ... (m-1), to produce a set of n, n an 
even positive integer, output audio data samples V[i], where i = 0, 1, ...(n-1), comprising 
20 the steps of: 

a) calculating an array of stim data S^^dM according to 

S^ddM = S[k] + S[m-l-k] for k = 0, 1, ...(mi/2-1) 

b) calculating an array of difference data SsusLk] according to 

SsubM = S[k] - S[m-l-k] for k = 0, 1, ...(m/2-1) 

25 c) calculating a first output audio data sample by a multiply-ac cumulate operation 
according to 

V[2i] = V[2i] + N[2i, k]*S^DD[k] for k = 0, 1, ... (m/2-1) 

where N[2i,k]= , J(32^20(2^.1)Tr| 
64 J 

d) calculating a second output audio data sample by a multiply-accumulate operation 
30 according to 

V[2i+1] = V[2i+1] + N[2i+1, k]*SsijB[k] fork = 0, 1, ... (m.'2-i) 
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where N[2i+1, k] = cos|-^^^^-I^^ 



l))i2k+l)Tz\ 



64 J 

e) and repeating steps c) and d) for i = 0, 1 , . . . (n/2-1) to obtain a ftill set of output data. 

The invention further provides a synthesis sub-band filter for use in decoding digital audio 
5 data, comprising a means for receiving or retrieving an input sequence of data elements 
comprising encoded digital audio data, a pre-calcuiation means for calculating an array of sum 
data and an array of difference data using selected data elements from the input sequence, and 
a transform calculation means for calculating a first sequence of decoded output values using 
said array of sum data and a second sequence of decoded output values using said array of 
10 difference data. 

The invention is described in greater detail hereinbelow, by way of example only, with 
reference to the accompanying drawings, in which: 

Figure 1 is a block diagram of major functional portions of an MPEG audio encoder; 
1 5 Figure 2 is a block diagram of major fimctional portions of an MPEG audio decoder; 

Figure 3 is a flow diagram of an MPEG decoding procedure; 

Figure 4 is a flow diagram showing a generalised form of a procedure according to 
the present invention; and 

Figure 5 is a flow diagram illustrating a preferred implementation of the invention. 

20 

Figure 1 is a block diagram illustrating the major components of an MPEG audio encoder 
circuit 2 t;onstructed in accordance with the aforementioned standards document. In- the 
figure, an input signal 4, comprising a pulse code modulated (PCM) signal having a 48 kHz 
sampling frequency and a sample size of 16 bits per sample, is provided as input to the single 

25 channel encoder 2. The input signal is first mapped from the time domain into the frequency 
domain by a sub-band filter bank 8. The resulting coefficients are normalized with scale 
factors which may be transmitted as side information. The coefficients thus obtained are then 
quantized and entropy encoded by a quantizer and encoding circuit 10. Masking thresholds 
of the quantization errors are calculated based on psychoacoustic values provided by a 

30 psychoacoustic model 14 to control the quantization step. The bit allocation is transmitted as 
side information. The coded signal is then multiplexed by a frame packing circuit 12 and an 
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encoded bitstream 6 is produced at the output of the encoder 2. 

A block diagram illustrating the main components of an MPEG audio decoder circuit 20 is 
shown in Figure 2. In the figure, an encoded bitstream 22 is provided to the input of the 

5 decoder. A bitstream unpacking and decoding circuit 26 performs an error correction 
operation if such operation was applied m the encoder. The bitstream data are unpacked to 
recover the various pieces of encoded information, and a reconstruction cncuit 28 reconstructs 
the quantized version of die set of mapped samples from the frames of input data. An mverse 
mapping circuit 30 transforms the mapped samples back into a uniform pulse code modulated 

10 (PCM) output signal 24 that reproduces the corresponding input signal which was provided 
to the encoder. 

The foregoing descriptions of the encoder and decoder are specific to the MPEG standard, and 
it is considered to be within^ the skill of those in the an to implement the various hardware 
15 functions described above. Accordingly, a more detailed hardware description of an MPEG 
coding system is not considered necessary for a full and complete understanding of the 
invention. It should be appreciated the invention described herein, although described in 
connection with the MPEG coding standard, is considered useful for other coding applications 
and standards. 

20 

Referring to Figure 3, there is shown a flow diagram 40 of steps involved in signal processing 
in layers 1 and 11 in an MPEGl audio decoder. To begin with, the bit allocation of an input 
bitstream (42, 44) is decoded (46). Thereafter, various scale factors are also decoded (48) and 
the samples are requantized (50). The encoded signal is decoded in a synthesis sub-band filter 
25 (52) and the decoded pulse code modulated signals are output (54, 56) for further processing 
and/or real time reproduction. The present invention relates primarily to the synthesis sub- 
band filter portion of the decoding process, when implemented for MPEG decoding. 

The synthesis sub-band filter bank is composed of two main functions, an Inverse Modified 
30 Discrete Cosine Transform (IMDCT) and an Inverse Pseudo-Quadrature Mirror Filter 
(IPQMF). The IMDCT, which can be viewed as an overlap transform, performs a 32 x 64 
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cosine mod-ulation txansformation, which means a frequency shift of a filter bank into one 
single filter. 

Consider a system in which output sub-band audio signal samples Vj ( i=0....63) are decoded 
5 from sequences of 32 encoded input samples S^, k = 0....31. The inverse MDCT of the 
sequence S^, is defined as follows: 




(1) 



for i=0,l, 63 



Taking the cosine symmetric property wJaerein: 



cos6=cos(2-n:-e) 



10 the IMDCT definition equation (1) may be modified as given below: 




(3) 



for f=0,l, 31 



This equation (3) may be computed according to the following algorithm: 



15 



repeat i = 32 times 

repeat k = 16 times 



if I is even, Sum = S[k] + S[31-k] 
if I is odd. Sum = S[k J- S[31-k] 
V[i] = V[i] + N[i, k] * Sum 



end k 



20 



end i 
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where 



k 



is the index of output samples (i = 0....3I) 
is the index of input samples (k=0...,15) 



N(i,k) = 




5 



S[k] 
V[i] 



represents the input sample data sequence 
represents the output oflMDCT 



The IMDCT equation, making use of tlie symmetrical property, is given in Equation (3) 
above, and the computational effort required for MPEG audio decoding is in large part 
dependant upon the efficiency wilh which the input samples can be processed through the 

10 IMDCT to obtain respective sub-band filter PCM samples. Embodiments of the present 
invention are able to reduce the number of arithmetic operations performed in implementing 
the IMDCT portion of the decoder, to thereby increase the computational efficiency of the 
decoding process. In particular, the number of addition operations required for the 
implementation of this equation can be reduced substantially by pre-computing the sum and 

15 difference of the sample data which is the input to the IMDCT. In addition, the pre- 
computation can take place outside the main IMDCT computational loop. Hence the main 
loop contains only the MAC operations, which can be executed very efficiently by any general 
purpose DSP in a minimum number of cycles. 

20 In the present invention, the dequantised sample data (e.g. 32 samples) from the encoded 
bitstream is pre-processed as per the symmetrical property of the cosine coefficients. The 
sample data is then split into two banks, each containing 16 samples. The sum and difference 
of respective data elements in the two banks is computed and stored in two arrays. These 
arrays are used as the input data for the subsequent MAC operations. 



Prior art implementations of equation (3) have required 32 x 16 Multiply-Accumulate 
operations and 32 x 16 Addition operations. By using the pre-computation operations 
described above, however, the number of Addition operations reduces to 2 x 16. This results 
in a saving of 30 x 16 Addition operations per Sub-band filter implementation, which in turn 
30 translates to a corresponding reduction in overall computational power. 



25 
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In the IMDCT equation (3), represents a sequence of m input data samples, where k = 
0 ... (m-l). In a typical implementation for MPEG decoding 32 input data samples may be 
processed, such that m=32. For pre-computing the sum and difference of respective data 
elements, the input data sample sequence is first arranged into two equally sized data banks, 
5 one constituting the high order data elements and the other the low order data elements: 
Data bank (1) for k = 0 ... (m/2)-l 

Data bank (2) Sk for k = (m/2) . . . (m-l) 

For example, in a preferred embodiment of the present invention where m=32, is split into 
10 two data banks comprising: 

(1) Sk for k = 0 .. 15 

(2) Sk fork = 16 .. 31 

The simi and difference data are calculated using respective data elements from the two data 
15 banks and is stored in two arrays of data, S^dd ^'^'^ ^svb. which are computed as follows: 
S^^^[k]=S[k]*S[m-l-k] _ for k = 0, I, (m/2) -I (4) 



Ssyg[k]=S[k]-Sim-l-k] for k = 0, I, (m/2)-l (5) 

In the aforementioned example of 32 input data samples, equations (4) and (5) reduce to: 
SaddM = S[k] + S[31-k] fork = 0, 1, ...15 

SsubM = S[k] - S[31-k] for k = 0, 1, ...15 

20 

The IMDCT equation (3) may now be divided into two portions and rewrinen as follows: 



T/r-i (32+/}(2A:+l)'n: „ 

64 "^^"^'^ (6) 
for /-0,2,4,....30 
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V[i] =X, cos-^ *55c;sW 



^^t^sL'^J (7) 
for /=1,3,5, 31 



As shown in the above equations (6) and (7), the IMDCT may now be calculated in two 
passes, an 'even pass' where the sum of the sample data is used (equation (6)), and an 'odd 
pass' where the difference of the sample data is used (equation (7)). The computational 
algorithms of the above equations are shown below. 

5 

Calculation of sum and difference vf sample data (Addition operations) 
repeat k = 16 times 

Sm^dK = + S,,^, 

10 endk 

Calculation of 'even' data of IMDCT (Multiply-Accumulate operations) 
repeat i — 16 times 

repeat k = 16 times 
15 Vfi] = VfiJ + N[i,k]*S^r^[kJ 

endk 

end i 

Calculation of 'odd' data of IMDCT (Multiply-Accumulate operations) 
20 repeat i = 16 times 

repeat k = 16 times 

Vfij = VfiJ + N[i,kJ*SsaBlk] 

end k 
end i 

25 

where i is the index of output samples (i = 0...3I) 

k is the index of input samples (k-0....15) 
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N(i,k) = 




5 



SfkJ 



VfiJ 



^ADD 



>SU3 



represents the input sample data sequence 
represents the sum of data arrar/ 
represents the difference of data array 
represents the output of the IMDCT 



Figures 4 and 5 illustrate the above procedure according to a preferred embodiment of the 
invention in the form of flow diagrams. The representation shown ia Figure 4, illustrates the 
general steps involved, and the procedure illustrated in the flow diagram 80 of Figure 4 

10 corresponds to the synthesis sub-band filter step 52 of the overall decoding procedure 40 of 
Figure 3. To begia with the input samples are received (82, 84) after having been isolated 
from the frames of encoded data received or retrieved. The input data samples are then 
utilised for pre-calculation of sum and difference data (step 86), as described above. This 
involves dividing the input data sample set into two equal sized sub-sets, which in the 

15 preferred embodiment consists of a fu-st sub-set comprising the lower order data and a second 
sub-set comprising the higher order data. For example, in the case of 32 input samples Sq to 
S32 as described, the first sub-set of input sample data may comprise the lower order input 
data So to 855 and the second sub-set comprises the upper order data samples to S31. 
Respective ones of each sub-set of input sample data are then used to obtain a sets of sum and 

20 difference data, S^dd and 85^3. As can be readily ascertained from the above description, in 
the preferred embodiment the calculation of the sum and difference data is performed using 
the lowest order samples fiom the first set with the correspondmg highest samples from tiie 
second set. For example, in the case of 32 input samples, the sum and difference data 
elements may be calculated as follows: 



25 



S^od[0] = S[0] -f S[31] 
S^dd[1] = S[l] 4- S[30] 
S,od[2] = S[2] + S[29] 



SsubEO] = S[0]-S[31] 
Ssub[11 = S[l] - S[30] 
Ssub[2] = S[2] - S[29] 



30 



Sadd[15] = S[i5] + S[16] 



Ssub[1-5] = S[I5] - S[16] 
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Once the arrays of sum and difference data have been calculated, the multiply-accumuiate 
operations required to calculate the IMDCT can be performed iteralively in two steps. The 
first step (88) is used to obtain half of the output samples (e.g. the "even" outputs) using the 
pre-calculated sum data comprising the S^dq data elements. The second step (90) is used to 
5 obtain the other half of the output samples (e.g. the "odd" outputs) using the pre-calculated 
difference data comprising the Sjua data elements. Each of these steps (88, 90) is an iterative 
multiply-accumuiate (MAC) operation involving each of the data elements from the respective 
S^DD or SsuB array. Furthermore, each of the MAC operations of steps 88, 90 are performed 
repeatedly (step 92) to obtain a fall complement of output samples. For example, where 32 
10 output samples Vg to V32 are required, each of the iterative MAC steps 88, 90 would be 
performed 16 times. Once the data for each output has been calculated, the data samples are 
output for PCM processing (step 94). 

A more detailed preferred embodiment of the decoding procedure is illustrated in the flow 
15 diagram 100 shown in Figure 5. Beginning at step 102, a sequence of m input samples S,, (k 
= 0 .... m-1) are received for decoding to n sub-band filter outputs V; (i = 0 .... n-1) at step 
104. In the preferred embodiment for an MPEG implementation, both the number of input 
samples m and the number of output samples n are the same, 32. Steps 106, 108 and 1 10 of 
procedure 100 for a loop for the pre-calculation process of determining and storing the sum 
20 and difference data arrays from the input data samples. The steps 112, 114, and 116 then 
form nested loops for the iterative multiple-accumulate calculation of the "even" ones of the 
output data^ elements (e.g. V; for i = 0, 2, 4, ... 14), using the pre-calculated sum data array 
Sadd- a calculation loop of steps 112 and 114 provides the iterative MAC operation, whilst 
the loop provided by step 116, enables calculation of each (even) alternate output data 
25 element. The remaining (odd) alternate output data elements are calculated in nested loop 
steps 118, 120, 122 using the difference data array S^^q. The resulting output sub-band data 
is then provided at fmal step 124. 

The preferred form of the invention presented herein results in a reduction of 480 addition 
30 operations per 32 sub-band samples. For a stereo output MPEG audio decoder, this is a 
reduction of 480 *36*2 arithmetic operations per frame. The overall reduction in arithmetic 
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operations which is achieved is approximately 46.875% per IMDCT. 

It will be readily apparent to those of ordinary skill in the relevant art that the present 
invention may be implemented in numerous different ways, without departing from the spirit 
5 and scope of the invention as described herein, and it is to be understood that such 
modifications are considered to be within the scope of the invention. In any event, it is 
immediately recognisable that one way the invention can be carried out, relating as it does to 
the processing of data, is using general purpose computing apparams operating under the 
instruction of software or the like which is produced separately and specially adapted to 

10 perform the methods of the invention. Alternatively, specialised computing apparatus such 
as a dedicated integrated circuit, chipset or the like may be constructed with the functions of 
the invention embedded therein. Many other variations to the particular implementation will 
of course be possible. It will also be recognised that in places in the description and appended 
claims where it is said that a data set is divided into sub-sets, for example, this division may 

15 be simply a notional one, and no physical separation need occur, as is known in the data 
processing art. 

The foregoing detailed description of the present invention has been presented by way of 
example only, and is not intended to be considered limiting to the invention which is defined 
20 in the claims appended hereto. 
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CLAIMS: 

1 . A method of decoding digital audio data, comprising the steps of obtaining an input 
sequence of data elements representing encoded audio samples, preprocessing the input 

5 sequence of data elements by calculating an array of sum data and an array of difference data 
using selected data elements from the input sequence, calculating a first sequence of output 
values using the array of sum data, calculating a second sequence of output values using the 
array of difference data, and forming decoded audio signals from the first and second 
sequences of output values. 

10 

2. A method as claimed in claim 1, wherein the array of sum data is obtained by adding 
together respective first and second data elements from the input sequence, the first and 
second data elements being selected from mumally exclusive sub-sequences of the input 
sequence. 

15 

3 . A method as claimed in claim 1 or 2, wherein the array of difference data is obtained 
by subtracting respective first data elements from corresponding second data elements of the 
input sequence, the first and second data elements being selected from mutually exclusive sub- 
sequences of the input sequence. 

20 

4. A method as claimed in claim 1, wherein the step of calculating an array of sum data 
and an array of difference data comprises dividing the input data sequence into first and 
second equal sized sub-sequences, the first sub-sequence comprising the high order data 
elements of the input sequence and the second sub-sequence comprising the low order data 

25 elements of the input sequence, calculating the array of sum data by adding together each 
respective data element of the first sub-sequence with a respective corresponding data element 
of the second sub-sequence, and calculating the array of difference data by subtracting each 
respective data element of the first sub-sequence from a respective corresponding data element 
of the second sub-sequence. 

30 

5. A method as claimed in claim 1, wherein the step of calculating a first sequence of 
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output values comprises performing a multiply-accumulate operation utilising each of ttie sum 
data elements. 

6. A method as claimed in claim 1 or 5, wherein the step of calculating a second 
5 sequence of output values comprises performing a multiply-accumulate operation utilising 

each of the difference data elements. 

7. A method as claimed in any preceding claim wherein the input sequence of data 
elements is derived from MPEG encoded audio data, and wherein the decoded audio signals 

10 comprise pulse code modulation samples. 

8. A method of decoding a sequence of m, m an even positive integer, input digital audio 
data samples S[k], where k = 0, 1, ... (m-1), to produce a set of n, n an even positive integer, 
output audio data samples V[i], where i = 0, 1, ...(n-l), comprising the steps of: 

1 5 a) calculating an array of sura data S^qdM according to 

SaddM = S[k] + S[m-l-k] for k = 0, 1, ...(m/2-1) 

b) calculating an array of difference data Ss^bM according to 

SsubM = S[k] - S[m-l-k] for k = 0, 1, ...(m/2-1) 

c) calculating a first output audio data sample by a multiply-accumulate operation 
20 according to 

V[2i] = V[2i] + N[i, k]*S;,DD[k] for k = 0, 1, ... (m/2-1) 

^ [(32+2i)(2^+l)ul 

where N[i, k] = cos -^^ —\ 

[ 64 j - 

d) calculating a second output audio data sample by a multiply-accumulate operation 

according to 

25 V[2i + 1] = V[2i+1] -i- Nli, k]*SsuBM ^ for k = 0, 1, ... (m/2-1) 

where N[i, k] = cos 



e) and repeating steps c) and d) 



(32H2/ + l))(2^^1)7t 
64 

"or i = 0, 1 , . . . (n/2-1) to obtain a full set of output data. 



9. A method as claimed in claim 8, wherein the number pi of input digital audio data 
30 samples is 32, and the number n of output audio data samples is 32. 
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10. A method as claimed in claim 8 or 9, wherein the decoding steps are repeated for 
decoding a series of frames of encoded audio data in an MPEG format. 

11. A synthesis sub-band filter for use in decoding digital audio data, comprising a means 
5 for receiving or retrieving an input sequence of data elements comprising encoded digital 

audio data, a pre-calculation means for calculating an array of sum data and an array of 
difference data using selected data elements from the input sequence, and a transform 
calculation means for calculating a first sequence of decoded output values using said array 
of sum data and a second sequence of decoded output values using said array of difference 
10 data. 

12. A synthesis sub-band filter as claimed in claim 11 wherein the pre-caiculation means 
and transform calculation means collectively perform an inverse modified discrete cosine 
transform of the encoded digital audio data. 

15 

13. An MPEG decoder including a synthesis sub-band filter as claimed in claim 11 or 12. 
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